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1 Introduction 



The calculation of the correlation functions of the spin chains and, in particular, their asymptotic 
analysis are very important problems in the field of quantum integrable models. In this article 
we consider one of the most representative examples of the lattice integrable models: the spin-^ 
Heisenberg chain. This model describes spin ^ particles situated in the sites of a one-dimensional 
lattice, interacting with their nearest neighbours, 

M 

H=Y. + ^m^m+l + - 1)) • (1-1) 

m=l 

Here A is the anisotropy parameter and dm^'^ are Pauli matrices, associated with each site of the 
chain. This Hamiltonian acts in a tensor product of M two-dimensional local quantum spaces 
Hm- ^ We imposed here the periodic boundary conditions 

For simplicity reason we consider here only the case when the number of sites M is even. 

Our main goal is to obtain some explicit expression for the two-point spin-spin correlation 
functions for the infinite chain at zero temperature. They can be defined as ground state mean 
values of products of local operators: 

5^^M =(^gK+ia^|^g), (1.2) 

a^-im) ={ipg\(y^+i(Ti\^g), (1.3) 

where | V'g ) is the ground state for the Hamiltonian (ll.lj) . 

To solve this model different techniques can be used. Originally the energy level were calcu- 
lated in 1958 by means of the coordinate Bethe Ansatz ^l^lElll]- Later (in 1979), an algebraic 
version of the Bethe Ansatz (or quantum inverse scattering method) was introduced by Faddeev, 
Sklyanin and Takhtajan Both methods permit to calculate the energy level and to obtain 
representations for the eigenstates and, in particular, for the ground state. However, for many 
years, computation of the correlation functions has been possible only for the free fermion point 

A = o 00121 cnj. 

First explicit results for more general situations, namely for the massive antiferromagnetic 
regime (A > 1), were obtained in 1992 by means of a completely different approach ((7- vertex 
operator method) by the Kyoto group jlH I12j. In 1996 a similar conjecture was written for the 
critical regime of the XXZ chain (— 1 < A < 1) ^31- The correlation functions were represented 
as multiple integrals with number of integrals equal to the distance. These results were confirmed 
(and generalised for the XXZ chain in a constant external magnetic field) in 1999 jl4l I15j by 

^AU the results of the present paper (with shght modifications) can also be applied to the XXZ model with 
external magnetic field. However for simplicity reason we consider here only the case of zero magnetic field. 
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means of the algebraic Bethe Ansatz and resolution of the quantum inverse problem |14[ I16j . 
This new approach permitted to understand better the results for the correlation functions and 
to obtain an asymptotic formula for a very particular correlation function (the so-called emptiness 
formation probability, which is the probability to find a ferromagnetic string of length m in the 
ground state) ^[E], and even an explicit result for this quantity for the point A = ^ [TO] . 

It is necessary to underline that the correlation functions calculated in 115j are not the 
two-point functions defined by (|1.2() and but they are the so-called "elementary blocks" 

(i.e. the mean values of products of local elementary 2x2 matrices with only one non-zero entry 
in consecutive sites). The two-point functions can be expressed as sums of these elementary 
blocks but the number of terms in this sum grows exponentially with respect to the distance 
(2'""'"^), which is very inconvenient for the computation of the long distance asymptotics. The 
main goal of the present paper is to reduce our expression for the two-point functions to a sum 
of m + 1 terms and thus to obtain a compact and manageable representation for it. One way to 
solve this problem was proposed in j2Uj . and here we present an alternative way to simplify the 
expressions for the spin-spin correlation functions. The advantage of this alternative approach is 
a simpler final form of the terms (each of them contains the same number of integrals and there 
is no need to introduce any additional integrations). It is also important to mention that each 
term in the final sum has a form very similar to the symmetric representation for the emptiness 
formation probability, obtained in [201, and thus there is a possibility to use the saddle point 
method similarly to JZj for its rough asymptotic analysis. 

The paper is organised as follows. In the next section we give a short reminder of the Algebraic 
Bethe Ansatz solution of the XXZ chain (following [S]) and an overview of our way to calculate 
the elementary blocks following our paper ■ In Section 3 we introduce the generating function 
of the correlation functions of the third components of spin and we show how to express this 
quantity in terms of the elementary blocks. We also explain how this function is related with the 
two-point function ()1.2|) . In Section 4 we prove our main result for the generating function and 
show how a similar technique can be applied directly to the correlation functions H1.3|) and ()1.2() . 
A further re-summation is performed in Section 5, which permits us to reduce the result to only 
one term. In the last section we show how this method can be applied in the free fermion point 
to obtain the known results for the two-point functions and their asymptotics. 

2 XXZ chain: algebraic Bethe Ansatz and elementary blocks 

First we have to determine the ground state of the Hamiltonian (jl.lj) . In the framework of 
the algebraic Bethe Ansatz it can be described in terms of the generalised creation and 
annihilation operators which can be obtained as elements of the quantum monodromy matrix. 
This matrix is completely defined by the i?-matrix of the model which, for the XXZ chain, is the 
usual trigonometric solution of the Yang-Baxter equation, acting in the tensor product of two 
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auxiliary spaces Fi (g) V2, Vi = C^: 



R{X) 



( 1 




V 





sinh A 
sinh(A— 
_ i sin ^ 
sinh(A— i^) 







i sin^ 

sinh(A— 
sinh A 
sinh(A— iC) 








1 / 



(2.1) 



the parameter C, being related to the anisotropy parameter^ 

A = cos C- 

The monodromy matrix can be constructed as a product of i?-matrices and can be written as a 
2x2 matrix in the auxiliary space: 



T{\) = Rom{X 



■C^ „ .X ^ .Cx / ^(A) B{X) 



(2.2) 



[0] 



Here {^} is a set of arbitrary inhomogeneity parameters and the i?-matrices Rok act in the tensor 
product of the auxiliary space Vq and the local quantum space Hk- The operator entries of the 
monodromy matrix A, B, C and D act in the same quantum space TC as the Hamiltonian of the 
XXZ chain. The commutation relations of these operators can be obtained from the Yang-Baxter 
equation: 

RuiX - fi) Ti{X) T2{fi) = T2{fi) Ti{X) RuiX - ^). (2.3) 

From this relation one can easily sec also that traces of the monodromy matrix taken in the 
auxiliary space (transfer matrices) commute for any values of the spectral parameters, 



[A{X)+D{X),A{ti) + D{fi)] = 0. 



(2.4) 



The Hamiltonian of the XXZ model can be reconstructed from the transfer matrix in the homo- 
geneous limit when all the inhomogeneity parameters are equal: = Oj 

d 



H = c-^logiA{X) + DiX)) 



+ const. 



(2.5) 



A=0 



It means in particular that the Hamiltonian commutes with the transfer matrix for any value of 
the spectral parameter, 

[H,A{X) + D{X)]=0, 

and the eigenstates of the transfer matrix (for arbitrary A) are eigenstates of the Hamiltonian. 

To construct the eigenstates of the transfer matrix one can use the operators B{X) as creation 
operator (and operators C (A) as annihilation operators) . This is possible if there exists a reference 

■^We give all the formulae in this paper for the critical regime (—1 < A < 1) of the XXZ chain, but it will be 
clear that similar computations can be done for the massive regime (A > 1) and for the XXX chain (A = 1). 
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state |0) which is an eigenstate of the operators A{X) and D{X) and is annihilated by the 
operators C(A) for any value of A. For the XXZ chain such a state exists and it is the ferromagnetic 
state with all the spins up. Now other eigenstates can be constructed by the action of operators 
B on this ferromagnetic state. More precisely, using the Yang-Baxter algebra 1)2. 3|) . one can show 
that the eigenstates of the transfer matrix (and hence of the Hamiltonian in the homogeneous 
case) can be constructed in the form: 

IV') = 5(Ai)...S(A^)|0), (2.6) 

where the spectral parameters {Aj} satisfy the Bethe equations, 

^ sinh(Aj -em + 4) TT sinh(Aj - Afc - zQ ^ . ^ , . 

iAsinh(A,-e™-z|)"L^-MA,-A, + <) ' ' ^-^^ 

The Bethe equations are very difficult to solve for a finite chain. However it can be shown that 
the ground state of the XXZ chain is described by this procedure in the homogeneous case with 
= and can be specified by a special choice of integers in the logarithmic form of the Bethe 
equations: 

M/2 

fc=i 

where the bare momentum po{X) and the "scattering phase" (?i>(A) are defined as 

sinh(A + 4) sinh(A + iQ 

Po{X) = ilog—— ^ (j){X) = i\og—— — . 

smh(A - 4) smh(A - iQ 

This state is the ground state of the XXZ chain in the homogeneous case, but even for the 
inhomogeneous model one can define a Bethe state with this choice of integers in the right hand 
side of the logarithmic Bethe equations. In the calculation of the correlation functions it will 
be convenient to consider such a state first and to take the homogeneous limit only in the final 
result. 

Even for the ground state there is no way in a generic situation to solve the Bethe equations 
explicitly. However, in the thermodynamic limit M — > oo, this state can be described in a very 
simple way in terms of the density of rapidities. More precisely: in the thermodynamic limit, for 
any smooth bounded function /(A), any sum over the Bethe roots corresponding to the ground 
state can be replaced by the following integral: 



M/2 t» 

^E/(^^)= / '^Ap(A)/(A) + o(^), (2.9) 

.7 = 1 ^ 
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where the density function p{X) can be obtained from a simple integral equation which replaces 
in the thermodynamic limit the Bethe equations, 



oo 

p(A)+ I dfip{fi)K{X-fi) = ^p'o{X), (2.10) 



where the kernel K{\) is a derivative of the "scattering phase" 



= ^"^'(^^ = ^ sinh(A+r)^sTnh(A-<) - ^'''^^ 
This equation can be easily solved by Fourier transform, 

2(; cosh(-^A) 

This information about the ground state is sufficient for the calculation of the correlation func- 
tions in the thermodynamic limit. 

The first problem which arises when one tries to calculate ground state mean values of prod- 
ucts of local operators in the framework of the algebraic Bethe Ansatz is the fact that creation 
(annihilation) operators B (C) are non local and do not permit a simple expansion in terms of 
spin operators. Thus it is rather difficult to establish commutation relations between these two 
types of objects. In our paper we proposed a way to solve this problem by expressing the 
local operators in terms of the monodromy matrix elements (hence solving the quantum inverse 
problem). The expression for the local elementary 2x2 matrices E^^"" with only one non-zero 
entry 

E'ji = Sje'Ske, 

can be written in a very simple form, 

E'y-= \{{A + D){iu-i\)T,^,.'JU-i\)\{{A + DY (2.13) 

k=l k=l 

It is easy to see that every local operator is expressed as a monodromy matrix element dressed 
with transfer matrices. In particular the operator cj^ is expressed as the dressed operator C{E,m — 
i|), (7~ as B{£^m — ^1) ^ ^(Cm " ^ 2) ~ ^i^m — '^f)- important to mention that this 

solution is very convenient to calculate the ground state mean values as the ground state is an 
eigenstate of the transfer matrix, 

{A{ik - ^h + D{i, - 4))| ) = n ^ c 

Now the correlation functions can be expressed only in terms of the monodromy matrix elements. 
The first type of object we can consider using this approach are the "elementary blocks" i.e. the 
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ground state mean values of products of the elementary local matrices in m consecutive sites, 

m 

F^{{eJ,e'^}) = {^l.g\llE'/'''\i;g). (2.14) 
i=i 

Using the solution of the quantum inverse problem such quantities can be expressed only in terms 
of the monodromy matrix elements, 

M/2 m M/2 

smh(Aj - ^fc + \ j=i k=i 3=1 



Fm{{e,,e'^])=\\{\{ 



^,;^,sinh(A.-e.-.|)y ,o,Wc(A/ff.(A.),0 



(2.15) 

(it is important to mention that operators B{\) are not normalised and thus to obtain the mean 
value, one should divide the r.h.s. by the norms of the Bethe vector corresponding to the ground 
state). Now we can use the Yang-Baxter algebra i.e. the commutation relations between the 
monodromy matrix elements. In particular we can act with the monodromy matrix elements on 
the dual Bethe state constructed by the actions of operators C(A). 

It is easy to see that after acting with all the operators in (|2.15|) on the dual Bethe state 
one obtains again a sum of states constructed by the action of the operators C(A) on the dual 
ferromagnetic state (but the spectral parameters do not satisfy any more the Bethe equations). 
Then, using the Gaudin-Korepin formula for the norm of the Bethe states [231 1241 I25j . and 
the determinant representation for the scalar products of a Bethe state with an arbitrary state 
|26| I14j. one can express the correlation functions (|2.14)) as sums of determinants. Now we have 
an explicit formula for the elementary blocks in term of the ground state solution of the Bethe 
equations. In the thermodynamic limit this sum can be simplified as all the sums over the Bethe 
roots can be replaced by integrals with density. The final result for any elementary block can be 
written as multiple integrals 

oo oo 

Fmiiej.e'j}) = T-r . TT / d\i... / (iAmJ'^({Afc},{^j,ej,e'})detm5'({Aj},{^fc}). 

j>k -oo -oo 

(2.16) 

Here the mxm matrix S does not depend on the choice of local operator and is defined uniquely 
by the ground state. Its elements can be written in terms of the density function (|2.12|) 

Sjk = p{Xj-^k). (2.17) 

The algebraic part J^({Xk}, {^j, ej, e'j}) arises from the commutation relation of the monodromy 
matrix elements and does not depend on the ground state. The expression for this function for 
the most general case can be found in the paper jJH]. Here for the calculation of the two-point 
functions we need only some particular blocks and we will give an explicit expression for the 
corresponding algebraic parts in the next sections. 
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Once the elementary blocks are calculated, any correlation function can be written as a sum 
of such multiple integrals. However the sums which appear for the two-point functions contain a 

number of terms that grows exponentially with the distance and in such a form cannot be used 
for the asymptotic analysis. For example, using the solution of the quantum inverse problem 
for the correlation functions of the third components of spin gzz{^)-> we obtain the following 
expression: 

smh A,-^, + . ^ ^^^^^^ _ . _ ^^^^ _ r 




It is easy to see that this function can be rewritten as a sum of 2'""'"^ elementary blocks. Similarly, 
for the function (m) , one has 

(m+lM/2 , 
smh(A,-a + zi 
y^^y sinh(A, -a-4) 

X { i^g mi - A) ( X{{A{ij - 4) + D{i^ - 4)) ] B{U+i - 4)1 ^9 )■ (2-19) 



and it can be written as a sum of 2"^~^ elementary blocks. The main goal of this paper is to obtain 
a manageable expression for these correlation functions from the multiple integral representation 
for the elementary blocks. 

3 Generating function 

In this section we consider the correlation function of the third components of spin ^^^(m). To 
calculate this function it is convenient to introduce a new object: the generating function 

m ^ 

Q^{l3)={^l;g\e^vWl,m]\^g), Ql,m = ^2^^"''i)- ^^"^^ 

To obtain the two-point function from the generating function one should take its second deriva- 
tive on (3 and second lattice derivative on m 



gzz{m) = (^2P^^ - 4P„— + 1 J 

where we used the standard definition of the first and second lattice derivatives 

T>mf{m) = f{m + 1) - /(m), = f{m + 1) - 2/(m) + f{m - 1). 



(3.2) 

/3=0 
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Using the solution of the quantum inverse problem we can rewrite this quantity in terms of the 
monodromy matrix elements: 

The generating function can be expressed as a sum of 2"* elementary blocks containing only 
diagonal elements of the monodromy matrix (operators — i|) and — i^)). For such 
elementary blocks the algebraic part can be written in a factorised form, as a product of "two- 
particle contributions" , 



smh(Aj - Cfe - 4) ' ^ ^ 



1 



oo oo 

J dXi... J dXmdetmSi{Xk},{^j}) 



n sinh(^j - ^fc) 

j<k -°o -oo 

TT sinh(Aj - gfc + iejO sinh(Afc - - i^kO .o .^v 

n sinh(A,-A, + % + efeK) ' ^'^ 

where numbers ej = ±| depend on the choice of the corresponding monodromy matrix elements 

It is clear that the elementary blocks with fixed number of operators of each type have a rather 
similar structure. Thus it is quite natural to put together the terms with the same number of 
operators D and to represent the generating function as power series on e^, 

m 

Qm{P) = Y.e'^F,im). (3.5) 

s=0 

The main result of this paper is a simple and compact formula for the terms Fs{m). This function 
can be expressed in terms of elementary blocks, 

\fe=l ,=1 Sinh(A,- - 6 - I a^^-^^-m=s 2 2 

(3.6) 

and hence the sums over all possible positions of operators D and A also can be written as 
multiple integrals 

oo oo 

Fsim) = \\ . [ dXi... [ dXmGs{m,{Xjmk})detmS{{Xj},{^k}), (3.7) 

fA smh(^j -ik) J J 

J^'^ — oo — oo 
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where the function Gs{m, {y^j}\{Ck}) under the integrals is a sum over the permutations a of the 
set 1, . . . , m 

r (rr. S\ ^^_ 1 ,^^-p. sinh(A^(j)-gfc + ie,(^.)C)sinh(A^(fc)-g,--ie,(,,)C) 

C.s(m,iA,||Ufc|j-^,^^_^^,^( 1) n sinh(A.(,)-A.(,)+z(e.(,)+6.(,))C) 

(3.8) 

where (— 1)^*^^ is the sign of the permutation a, the factorials in the denominator arise from 
the additional summations over permutations of the variables of the same type ( with ej = ^ or 
ej = — We set for j < s, ej = and for j > s, ej = ^. We use here a very important property 
of the multiple integral representation, namely, the fact that the determinant of densities does 
not depend on the choice of local operators. 

The function Gs{m, {Xj}\{^i^}) is a rational function of e'^^ and e^*^. It is by definition skew- 
symmetric under the permutations of Ai . . . , Am- The first property of this function which we 
need for our computation is its symmetry on the variables ^j. 

Lemma 3.1. The function Gsim, {Xj}\{^k}) defined by iS. <Sj) is symmetric under the permuta- 
tions of the variables (,i,^2, ■ ■ ■ ,£,m- 

Proof: To prove this lemma it is sufficient to show for any k that 

Gsim,{Xj}\^i,. . . ,^k,(,k+i,- ■ • ,^m) - Gs{m,{Xj}\(,i, . . .,(,k+i,Ck, ...,Cm) = 

Consider first this difference for corresponding monomials in 1)3. 8|) and note that two terms differ 
only in one "two-particle contribution": 

-TT sinh(Aj - 6 + ^ejC) sinh(A; - Cj - ieiQ yr sinh(Aj - + i^jO sinh(Afc - - iCfcC) 

j-^-L sinh(Aj - A/ + i{ej + e;)C) .^^ sinh(Aj - A^ + i{ej + ek)C) 

i^k 

sinh(Afc+i -ik + i^k+iC) sinh(Afc - jk+i - if^kC) 

sinh(Afc+i - Afe + i{ek+i + ek)C) 

sinh(Afc+i - Ck+i + ^efc+iC) sinh(Afc - (k - i^kC) 

sinh(Afc+i - Afc + i{ek+i + ek)C) 

t ^ TT sinh(Aj -^1+ iejC) sinh(A; - - ieiQ 

■■ smh(^fc - ^k+l) I I r-TTT r — —, ■ rr: 

-"-J- smh(Aj - Xi + i{ej + e/)C) 



l^k 



m 



^ , , sinh(Aj - gfc + iejC) sinh(Afc - - iekQ 



n 



j-=fc+2 sinh(Aj - Xk + i{ej + ek)0 

Consider now the same term for the following permutation of the set Ai, . . . Xyu- 

{Ai, . . . , Afe_i, Afc+i, Xk, Xk+2, • • • 1 Ami- 
It is easy to see that for this permutation we obtain exactly the same contribution with an 
opposite sign (as the sign of this permutation is —1) and the sum of these two contributions is 
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zero. Now we should take a sum of such monomials over all possible permutations of the set 
Ai,...,Am, but as this sum can be split into such pairs with permuted Xcr{k) and Xa{k+i) we 
immediately obtain that this sum is zero. □ 

Using this symmetry we can obtain recursion relations for the function Gs(m, {Aj}|{^fc}) in 
the points Xj = S,k — more convenient now to extract the common denominator and to 

consider the function {Aj}|{^fc}): 



r( r \ iirg i \ _ 1 TT sinh(Aj - Afc) 

s[m, 1 _ _Q ^.^^^^^ _ ^ ^^^^ ^ ^^^^^ ^.^^^^^ _ _ ^^^^ ^ ^^^^^ 

xG,(m,{A,}|{efc}) 



Directly from this definition we can establish two evident lemmas: 

Lemma 6.2. The function e^""-^^^^ Gs{m, {\j}\{Ck}) is a polynomial function of e^^^ of degree 
m — 1. 

Lemma 3.3. 

Go(l,Ai|ei)=G'i(l,Ai|ei) = l. 

These lemmas mean that if we obtain recurrence relations for this function in the points 
Xj = Cfc ~ ^2 then they will define it completely (as it is a polynomial of degree m — 1 defined in 
m points). The recursion relation for this function can be written in the following form 

Lemma 3.4. Gs{m, {Xj}\{(,k}) satisfies the following recursion relations, 
Gsim,{Ximk}) 



m 



sinh(Aj - - Jl sinh(Aa - 6 - 



X Gsijn — 1, Ai, . . . , Aj_i, Aj_|_i, . . . , Am|^i, . . . , ik~\iik^\^ ■ ■ ■ Cm), = (3-9) 



2 



Gsim,{Ximk}) 



Y[ sinh(Aj -Ca- Jl sinh(Aa - Ck - i-) 



X Gs-i{m — 1, Ai, . . . , Aj_i, Aj+i, . . . , Am|^i, . . . , Cfc~ii Cfc+ii ■ ■ - Cv 



1 

~2 

(3.10) 



Proof: Prom the Lemma l3. II it is clear that it is sufficient to prove these recursion relations 
only for Xj = — ^§ for ej = ^ and for Xj = Cm — ~ ~\- ^^^^ *° t^isX only 

terms with a{j) = 1 {cy{j) = m ) will survive in the sum over permutations ()3.8p . We obtain 
a sum over permutations of the set {Ai, . . . , Aj_i, Aj+i, . . . , Am} and using the definition of the 
function Gs{jn,{Xi}\{S,k}) we obtain immediately the recursion relations 1)3. 9p . 1)3. lOp . □ 
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These recursion relations define completely the function Gs{'m,{Xi}\{S^k})- It means that if 
we can find a function satisfying Lemmas 3.1-3.4 it is the function Gs{'m,{Xi}\{(,k}) we need. 
Such properties (without dependence on s) were established for the first time by Korepin in [2SI 
for the partition function of the six vertex model with domain wall boundary conditions. In fact 
the recursion relations obtained here are exactly the same as the Korepin's ones. The (unique) 
solution for these relations, found by Izergin in (which does not depend on s), satisfies the 
Lemmas 3.1-3.4 for any value of s and thus gives the function Gs(?7i, {A;}|{^fc}): 

Theorem 3.1 (Izergin 87). The only function satisfying lemmas 3.1-3.4 is 

where Z,m({A}|{^}) is the partition function of the inhomogeneous six vertex model with domain 
wall boundary conditions, which can he written in the following form, 

mm 

n n sinh(Aj - 6 + 4) sinh(Aj - Cfc - 4) 

Zm{{Ximu}) = ''^''^' ■ ^^■^.<c 7\ det™-M({A/}|fe}), (3.11) 

11 smh(Aj - Xk) smh(^fc - ij) 
j>k 



sin C 

sinh(Aj -^k + 4) sinh(Aj - Cfc - 4) 



Mjk = 7^ ^. (3.12) 



Proof: To prove this theorem one needs just to verify that the function Z„i({A}|{^}) 
satisfies all the Lemmas 3.1-3.4 for any value of s. Then, as it is a polynomial of degree m — 1 
coinciding with another polynomial of the same degree in m points, one can conclude that these 
two functions are equal. □ 

The fact that the function Gs{m, {A/}|{^fc}) does not depend on s is a very important pecu- 
liarity of this re-summation technique. Hence, the generating function can be written as a sum 
of m + 1 multiple integrals of the following form, 

oo oo 

^•("'> Vm^^C n imh(4F&) I"'^-^ I f«» 



-oo 



xe^(Ai,...,A^) Z^({A}|{e}), (3.13) 
where the only factor depending on s which we denote 6m({A}) 

pvwx A ^ = 17 fr ^ rr sinh(Aj - Afc) 



X 



TT '^"^^^^ - (3.14) 



s>j>k> 



^ sinh(Aj - Afc + <) sinh(Aj - Afc - iC) 
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will appear in all our results. 

Thus we reduced the number of terms from exponential to polynomial order. This represen- 
tation is also interesting because of its unexpected relation with the partition function of the 
corresponding inhomogeneous six-vertex model with domain wall boundary conditions. We also 
suppose that this representation can be convenient for the asymptotic analysis of the two-point 
function gzz{m). 

It is also very important to mention that for this re-summation we manipulated only the 
algebraic part of the expression for the elementary blocks and hence it can be done in a similar 
way for the XXZ spin chain in a magnetic field. The result in this case has a slightly more 
complicated deteminant of densities and different integration contours but the algebraic part is 
the same and contains the partition function Zm({A}|{.^}). 

A rough asymptotic analysis of the generating function can be performed using a modification 
of the saddle point technique introduced in [2111 for the emptiness formation probability. It shows 
that there is no gaussian contribution to this correlation function and that the main order can be 
written as C exp(^), as it should be, but it is not sufficient to describe the power-like behaviour 
of the two-point function. The interesting pecularity of this saddle point analysis is the fact that 
the "saddle point density" here coincides with the ground state density p{X). 

4 Two-point functions 

A similar re-summation can be done for the two-point functions g-\ (m) (or for gxx{m) and 

gyy{m)) and gzz{m)- The function g^ (m) can be written in the following form in terms of the 

monodromy matrix elements: 



To calculate this function one should sum up 2*" elementary blocks of the following form 
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which can be written as multiple integrals in a factorised form 

-TT TT- smh(Aj - Cfc + 
,y y sinh(A, -a-4) 



m+1 

n sinh(^j -^fe)-oo 

k>j>l 



J d\2... j dXm J dX+ J dA_det^+iS'({A2,...,A^,A+,A_},{^j}) 



sinh(Aj - Cfc + iejO sinh(Afc - - ie^C) J sinh(A+ - 6 + 4) sinh(A_ - - i^) 



X 



•>j^.>2 sinh(Aj - Afc + i(ej + efc)C) / sinh(A+ - A_) 

— sinh(A_ - 6 - 4) sinh(Afc - 6 + 



n 



sinh(A_ - Afe + z(ejk - i)C) 



-TT sinh(A+ - Cfc + 4) smh(Afc - C^+i - iekQ 
l}^ sinh(A+ - Afe + i{ek + ^)C) ' ^ 



The blocks with the same number of operators D (ej = can be put together 

m— 1 

g+-{m) = g+-im, s), (4.3) 



s=0 

where 

02 H ham— m+l=s 

It is easy to see that in all the terms of this sum written as multiple integrals the determinant of 
densities and all the factors containing A+ and A_ arc the same. Thus to obtain the corresponding 
algebraic part it is sufficient to take the sum over all possible permutations of A2, . . . , A,„ of the 
product 

sinh(Aj - ^fc + i€jC) sinh(Afc - - if^kC) 



n 



j>k>2 sinh(Aj - Afc + i{ej + ek)0 

This sum has exactly the same form as the corresponding term for the generating function. 
Using exactly the same arguments as in the previous section we obtain for the contributions to 
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the two-point function: 

oo oo oo oo 

5_| (m,s) = — ^ — IT . , ,^ / dXo ■ ■ ■ / dXm / dA+ / dX- 

sinh(A_ - 6 - 4) sinh(Ajk - 6 - 4) 1 



n 



yk=2 sinh(A_ - Afe - iQ J 

' Y\ sinh(A- - ^fc - 4) sinh(Afc - + i|) 

^TT^ sinh(A+ - Cfc + 4) sinh(Afc - + i^) 
^li sinh(A+ - Afe) 



/ 1^ sinh(A+ - ^fc + i|) sinh(Afc - jm^^ - i|) \ 
Vii sinh(A^-A. + ^C) ) 

sinh(A+ - gi + i|) sinh(A- - - i\) n x ^ 

sinh(A+ - A_) •B^_i(A2,...,A^) 

xZm-l{{^2, ■■■ , Am}|{^2, • • • , ^m}) det^rj+l <S'({A2, . . . , A^, A+, A_}, {^i, . . . , ^^+1}). 

(4.5) 

This result is shghtly more complicated but very similar to the result for the generating 
function. 

One should note that a very similar formula can be written for the correlation function 
gzzi'm) directly (without any use of the generating function). This function can be written in 
the following form in terms of the monodromy matrix elements: 

m+lM/2 . , 

smh(Aj - ^fc + 



5..(m)=-i+4 n n 



k=i j=i sinh(Aj - ^fc - q) 
To calculate this function one should sum up 2'"^^ elementary blocks of the following form, 

(m+lM/2 . , 
k=i ,=1 sinh(Aj- - Cfe - i^) 



U=2 
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As usual we put together the blocks with the same number of operators D, 

m—l 



gzz{m) = + gDDim,s), (4.7) 



s=0 

where 

m+lM/2 , /- , -Cx 

smh(Aj - U + ^2) 



9DD{m,s) = n n • . . -C 

' smh(Aj - 6 - 4 



a2H hflm— m+l=s 



(4.8) 

Using exactly the same arguments as for the function (m) we obtain for the contributions 

to the two-point function: 



oo oo oo oo 



gDD{m,s) =— . TT ■ , TT / dX2... dX^ / dX^+i 

s!(m - 1 - s)!sm™' V J- , ^ smh(4,- - i II II 

^ ' m+l>j>fe>l ■' 



-oo — oo — oo — oo 



/ sinh(Ai - 4 + 4) sinh(Afc - gi - 4) \ 
""[l}^ sinh(Afc-Ai-<) J 

/ sinh(Ai - gfc + 4) sinh(Afc - 6 + 4) ^ 

( sinh(A,n.+i - gfc - 4) sinh(Afc - jm+i + 4) \ 
sinh(Am+i - Afe - iQ j 

I sinh(Am+i - gfc - 4) sinh(Afc - ^m+i - 4) ) 

[k=s+2 sinh(A^+i-Afc) J 

sinh(Am+i - - 4) sinh(Ai - Cm+i + 4) r>s / x x x 

smh(Am+i - Ai - <) 

X^m-l({A2, • • • , Xm}\{C2, ■ ■ ■ , Cm}) det^+l S'({Ai, A2, . . . , Am, X-m+i}, {,^1 , . . . , Cm+l})- 

(4.9) 



5 Further re-summations 

In this section we show how one can proceed to a complete re-summation of terms and reduce 
the result for the two-point correlation functions to only one term (written again as a multiple 
integral). To do it one should note that all the terms in the sum p.5() have a very similar struc- 
ture. Extracting the common denominator one can obtain the following sum for the generating 
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function: 

oo oo 



^ s\{m — s)\ sm C sinn(4j — i ' ' 



" ^ / / 

oo — oo 

mm 

n n sinh(Aj - Cfc + 4) sinh(Aj - ^fc - i§) 
fc=i j=i 

^ i 

W sinh(Aj — Afc) sinh(Aj — A^ + iC,) sinh(Aj — A^ — iQ) 

j>k 

X ZUWm) det^S{{X},{^})Hs{{X}m)Hs{{X}m), (5.1) 



where Hs means complex conjugation and the function ff^dA}!!^}) is defined as 

s m 

n n sinh(Aj - Afc + <) H sinh(Aj - A^) sinh(Aj - A^) 

^ r \ 1 I rr-l \ j=lfc=s+l l<i<'c<s s+l<j<A:<m 

-'^Hi^lli?!) = 7 ^ 

mis m \ 

n n sinh(Aj - Cfc - if) n sinh(Aj - 6 + i§) 



(5.2) 



It is easy to see that this function can be written as a Cauchy determinant 

As) 



g.(Wltt})=n ,,.i.fe_;^) ''°'« 

sinh(Aj — — 4' 



As) 



1 



^ sinh(Aj - Cfc + 4) ' 

The sum over s in ()5.1|) can be taken under the integrals. Here we separated the terms which 
depend on s: 

m jjg 

V —-^ det det W(^). 

^-^ s'.im — sy. 

This sum can be simplified if one introduce some auxiliary contour integrals: 

/3s 

y — ^ TT det W(^) det W(^) 

s=o ^ ^ 

1 / dzi / / j3 I / TT sinh2z 



— I r . . . ^ ■ exp I . / V -7 I " ; I I I I A A 

m\ J 2z7r / 2«7r \ *C 2 / I sinh(zj - 4) sinh(zj + 

where 



n . C L T-A det^+det^-, 



J"; 



1 



sinh(Aj ± - ^k) ' 
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and contours are chosen in such a way that the points inside the contours and 

all the poles at the points Xj it — ^fc = are outside. Thus we can represent the generating 
function as a single term but the number of integrals is 2m now: 



oo oo 

dz^ f dzr^ 



1^'' — oo — oo 



m m 



n n sinh(Aj -Ck + 4) sinh(Aj - - 
k=lj=l 

Y\ sinh(Aj — Afc) sinh(Aj — A^ + iC) sinh(Aj — A^ — iQ 

j>k 



det„S({A},U}) 



™ exp \^^{zj + sinh22:j 
sinh(2;j — «|) sinh(2;j + z| 



n ^ ^''^ ^"^^ ^~ ■ (^•^) 



Now, using the symmetry of the expression under the integral with respect to permutations of 
pairs {Xj, zj), we can replace one of the determinants det^^ by a product of its diagonal terms. 
It permits in particular to separate the variables zj and hence to integrate over them. It is easy 
to see that it leads to the following final result for the generating function: 



oo 



QM=-^ll . , , [dXi...[ dA™Z„({A}|{e})det„5({A},{e}) 



j<k 

m m 



n n sinh(Aj - Cfc + 4) sinh(Aj - Ck - 4' 

detmG, (5.4) 



Y\ sinh(Aj — Afc) sinh(Aj — A^ + iQ sinh(Aj — A^ — iC) 

j>k 

where the m x m matrix G is defined as 

1 

Qjk = 7 T- H 7 (5.5) 

sinh(Aj — ^k + 4) sinh(Aj — ^j — i |) sinh(Aj — Ck — 4) sinh(Aj — ^j + z|) 

Similar results can be obtained for the two-point functions. This is the most compact formula 
for the generating function (we reduced the number of terms from 2™ to one). However we think 
that the formulae obtained in the two previous sections are more convenient for the asymptotic 
analysis as they permit a natural homogeneous limit (which is not the case of 1)5. 4() ). 



6 Free fermion point 

As the first application and check of this new re-summation formula we consider the free fermion 
point (C = f )• Of course the representations for the correlation functions in this point have 
already been obtained by different methods, but the formulae and 1)4.51) give a simple and 

elegant way to get these explicit results. 
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6.1 Generating function 

We calculate the generating function 



m ^ 



This function can be written as 



s=0 



(6.1) 



where contributions Fs{m) are given by (|H.1I-{|) . 

Taking into account that ^ = ^ all the determinants can be calculated. In the homogeneous 
limit one obtains the following multiple integral representation: 



FJm) 



7r™'s!(m — s)! 



dXi . . . dXm sinh^(Aj — 

Joo l<j<k<s 



X Y[ sinh2(Aj - Afc) Jl Jl cosh2(Aj - A^) JJ cosh"™ 2Xj 

s<j<k<m j=l k=s+l j=l 



(6.2) 



After changing variables 



= -log(tanpj), 



we obtain a much simpler representation: 



FJm) 



vr^-sllm — SI! 



2 2 

dpi... I dpm Y\ sin^i^ jPj - £kPk)- 
j>k 



(6.3) 





The expression under the integral can be rewritten as a product of two Vandermonde determi- 
nants: 

2 



Y[ s'm^ [ejpj - EkPk) =2" 
j>k 



j>k 



2"^-™ detV{{£p})detV*{{ep}), 



where star means hermitian conjugation. The product of these two determinants can be calculated 
as the determinant of the product of two matrices: 

detV{{ep})detV*{{ep}) = det{V* {{sp})V {{ep})) 

m 

{V*V)kn = ^e'^'^^-'''^'^PK 
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Now we should take into account that we integrate this determinant. It can be written as a sum 
over permutations, 

2 2 

-^^("^) = ^TTT TT / dpi-- - / dp.myidetVF(cr,s,{p}), 

7r"^s\[m — s)\ J J ^-^ 



where the matrix W is defined as 

As we integrate over all the variables p, permutations of the variables with the same value of 
the parameter e do not change the integration result and thus the sum over permutations can 
be replaced by a sum over the partitions of the set {p} into two subsets {p"*"} and {p^}, with a 
number of elements in the first one being s, 

ZL IL 
2 2 

= [ dpi-- - [ dprn J2 W{{p^}, {p-}). 

Here the matrix {p~}) is defined as 

e„ being + or -. 

Now considering the entire sum over s (|6.1j) . one can easily note that it can be rewritten as 
a determinant of a sum of two matrices: 

— — 

2 2 

2m (/^) j dpi--- j dpm detm U (P, {p}), 



Ukn{f3, M) =e2^('=-")P" + e/3e2^("-'=K . (6.4) 
It is now possible to calculate all the integrals and to write the final result as a determinant: 
QmiP) =detmT{p), 

TUP) =Skn^- + (1 - 4n)(l - e^) 2iX-k) • ^^-^^ 
From this formula one can easily obtain the well known result for the two-point function gzz {fn) : 

Ka^+i) = ^((-ir-l). (6.6) 
The emptiness formation probability 29 comes directly from (|6.5|) by taking the limit /? — > — oo. 
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6.2 Two-point functions 

The two-point functions {(^i^'l^j^i) and (cfc^+i) '^^^ calculated in a similar way in the free 
fermion point. 

First of all a re-summation formula can be obtained for these functions almost in the same 
way as for the generating function. We can consider even a more general function: 

9+-{m,(3) = {af exp(/?Q2,m)o-~+i), 

m 

where Q2,m = i(l ~'^|)) which includes also the correlation function of fermionic fields. After 
re-summation it can be represented as: 

m—l 

g^.{m,P) = Y,e'f'9t-{rn), (6.7) 

s=0 

where the contributions (m) are given by (|4.5j) . This general re-summation formula can be 
simplified in the free- fermion case: 

„ oo oo oo oo 

— oo — oo — oo — oo 

s+1 m 

n sinh2(A, -Afc) n sinh2(A, -Afc) n FI cosh2(A, - A^) 

2<j<k<s+l s+l<j<k<m j=l k=s+2 

n cosh" 2Xj 

s+1 m 

cosh(A+ — A_) n sinh(A_|_ — Xj) cosh(A_ — Aj) W cosh(A+ — Xj) sinh(A_ — Xj) 

i=2 j=s+2 

^ cosh'"+i 2A+ cosh^+i 2A_ 

X sinh'"(A+ - i^) sinh'"(A^ + i^). (6.8) 

Changing variables as in the previous case 

= ^log(tanpj), 

we obtain: 

TT TT TT TT 

9 , "2 "2 "2 "2 

— 1 \s r f f f 

9t~irn) = .^_,,^.^,„„^ dp,... dp^ dp^ / dp^e^(^+-^-)- 



z™-i7r™+is!(m - 1 - s)! 





m 

X Jl s\v?{ejPj - EkPk) JJ sin(p+ - SjPj) sin(p_ + SjPj) sin(p+ + pJ). (6.9) 

m>j>k>l j=2 
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This expression also can be rewritten as a product of two Vandermonde determinants, but now 
the matrices have different sizes: 

"2 "2 "2 "2 



X detm+1 V{m + l,{p+, -p-,e2P2, ■ ■ ■ ,emPm}) detm-i V*{m - l{e2P2, • • • ,emPm}), 



The product of these two determinants can be again rewritten as a determinant of a product of 
two matrices if we add two rows and two columns to the second one: 

mm + 1, M) =V*,{m + 1, M) = 6i„ 

Vijim + 1, M) =V*2{m + 1, M) = <52„ 

VCjim + 1, M) =V,Uj-2im - 1, M), A:, i > 2. 

Now we can take the product of these two determinants and proceed with the sum over s in the 
same way as for the generating function: 

TT TT TT TT 

2 2 2 2 



f/fcn(/3, W) =e2^('=-"-i)P" - e/3e2^("-fc+iK. (6.10) 
Now all the integrals can be calculated. For the two-point function {(3 = 0) we get 

fi'+-M = ^m+i det^+ir~+(m), 

rp-+ _rp-+ _ 

-'ii —^21 — 2 ' 
1" 2„ 2(n-l) ' 



Tkn" =T^^^Tr^- (6-11) 



1 + (-1)^-" 
(k-n-l) 

This determinant can be easily computed and gives a formula obtained by Wu in 7^ (a more 
general formula was obtained by McCoy 0): 



[m + l] 



+ _ ,_ {-l)"' -fi T^jk) 'fj' T^jk) 
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The asymptotic behaviour of such a product was also obtained in [3 using the technique intro- 
duced in [HO] : 



+ 
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